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Abstract. We report on total-energy electronic structure calculations in the density- 
functional theory performed for the ultra-thin atomic layers of Si on Ag(lll) surfaces. 
We find several distinct stable silicene structures: x ^/3, 3x3, ^/7 X \/7 with the 

thickness of Si increasing from monolayer to quad-layer. The structural bistability 
and tristability of the multilayer silicene structures on Ag surfaces are obtained, 
where the calculated transition barriers infer the occurrence of the flip-flop motion 
at low temperature. The calculated STM images agree well with the experimental 
observations. We also find the stable existence of 2 x 1 7r-bonded chain and 7x7 dimer- 
adatom-stacking fault Si(lll)-surface structures on Ag(lll), which clearly shows the 
crossover of silicene-silicon structures for the multilayer Si on Ag surfaces. We further 
find the absence of the Dirac states for multilayer silicene on Ag(lll) due to the 
covalent interactions of silicene-Ag interface and Si-Si interlayer. Instead, we find a 
new state near Fermi level composed of tt orbitals locating on the surface layer of 
\/3 X \/3 multilayer silicene, which satisfies the hexagonal symmetry and exhibits the 
linear energy dispersion. By examining the electronic properties of 2 x 1 7r-bonded 
chain structures, we find that the surface-related tt states of multilayer Si structures 
are robust on Ag surfaces. 
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1. Introduction 

A honeycomb-structured two-dimensional (2D) atomic layer consisting of group IV 
atoms exhibits peculiar electronic properties due to a fact that electrons near the Fermi 
level {Ey) follow effectively the massless Dirac equation (Weyl equation) [1], A well- 
known and only unequivocally measured example is graphene where intriguing properties 
such as an anomalous quantum Hall effect is observed [2, 3]. Another element in group 
IV, Si, which has sustained our modern life, should exhibit such fascinating properties 
[4, 5] thus opening a new door to the next-generation technology with its pronounced 
relativistic effects related to the spin degrees of freedom [6, 7, 8]. 

The 2D structure of Si, known as silieene, is first synthesized on Ag(lll) surfaces. 
The monolayer silieene on Ag(lll) exhibits variety of structures. The superperiodicities 
of 4 X 4 [9, 10, 11], \/l3 X \/l3 [10, 11], and 2\/3 x 2\/3 [11, 12] with respect to 1 x 1 
Ag(lll) surface have been observed, and the simulated scanning tunneling microscopy 
(STM) images of theoretically determined structures reproduce the observed images 
excellently [13, 14]. Recently, the \/3 x silieene on Ag(lll) with respect to 1 x 1 
silieene has been reported by Chen et a/.[15, 16]. They also claimed the flip-flop 
motion between the two stable hexagonal \/3 structures to explain their honeycomb 
symmetry STM images at higher temperature and its freeze at lower temperature [16]. 
However, our density-functional theory (DFT) calculations have clarified that Chen’s 
observation should correspond to the bilayer silieene [17]. On the other hand, the 
multilayer (including bilayer) silieene structures have been also synthesized on Ag(lll) 
surfaces in other experiments [18, 19, 20, 21, 22], which interestingly exhibit the similar 
honeycomb x STM images as that reported by Chen et a/., irrespective of the 
thickness of silieene. Thus it is very important to clarify the structural details for the 
multilayer silieene on Ag(lll). 

On the other hand, Si surfaces are premier stages on which the fabrication of 
most electron devices is achieved. Structural characteristics of such Si surfaces are well 
identified experimentally and theoretically [23]. In particular, the atomic layers of the 
(111) surface has the hexagonal network in the lateral plane, and thus the similarity and 
dissimilarity between the silieene and the surface atomic layers of Si(lll) surface are 
intriguing. Actually the cleaved Si(lll) surface shows 2x1 superperiodicity and then 
its annealing converts it to the 7x7 periodicity. The 2x1 structure has been identified 
as the TT-bonded chain structure [24] in which the top two Si layers are drastically 
reconstructed to form chains of tt orbitals associated with the five- and seven-membered 
rings of the first and the second surface layers. The 7x7 structure is more complicated 
and a dimer-adatom-stacking-fault (DAS) model [25, 26] is now established [27, 28]. It 
is important and interesting to clarify how a crossover between the silieene structure 
and the reconstructed surface structure takes place in the multilayer silieene. 

As for the electron states near Flp, the situation is controversial for the 
monolayer silieene on Ag(lll): Angular-resolved photoelectron spectroscopy (ARPES) 
measurements [9, 29] show the existence of the electron state with the linear energy 
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dispersion (Dirac state) for the 3x3 silieene on 4 x 4 Ag(lll), whereas the DFT 
calculations clarify the absence of Dirac electrons due to the strong silicene-substrate 
interactions [13, 14, 30, 31, 32]; it is of note that no Landau-level sequences peculiar 
to Dirac electrons have been observed [33]. From the interference patterns obtained by 
scanning tunneling spectroscopy (STS) measurements, Chen et al. claimed the presence 
of Dirac states with the linear dispersion [15, 16] for the \/3 x silieene on Ag(lll), 
whereas Arafune et al. deduced the contrary conclusion, the absence of Dirac electrons, 
from essentially identical STS experiments [34], Our DFT calculations show that the 
\/3 X silieene structure observed by Chen et al. actually correspond to the bilayer 
silieene and the linear energy dispersion is attributed to the sp band of Ag [17]. Recently, 
the ARPES measurements also observed the linear energy dispersions for the \/3 x '/S 
multilayer silieene on Ag(lll), where the electron velocity is much smaller than that 
observed by Chen et al. The linear energy bands are recognized to be the Dirac tt 
and anti-bonding vr (labeled as vr*) bands of multilayer silieene [19, 35]. Details of the 
electronic structure of the multilayer silieene are still mysterious and further theoretical 
work is in great need for the clarification. 

In this work, we report on total-energy electronic structure calculations for the 
ultra-thin atomic layers of Si on Ag(lll) surfaces. We find several distinct stable silieene 
structures: ^/3 x \/3, 3x3, \/7 x ^/7 with the thickness of Si increasing from monolayer 
to quad-layer. We also find the existence of 2 x 1 and 7x7 Si(lll)-surface structures 
on Ag(lll). The crossover of silieene-silicon structures on Ag surface appears when the 
thickness of Si layer becomes larger than 2. We also find the structural bistability and 
tristability of the multilayer silieene structures, where the calculated transition barriers 
infer the existence of flip-flop motion at low temperature. We further find the absence 
of the Dirac states for multilayer silieene on Ag(lll) due to the covalent interactions 
of silieene-Ag interface and Si-Si interlayer. Interestingly, we find a new state near E-p 
composed of tt orbitals locating the top layer of \/3 x ^/3 multilayer silieene, which 
satisfies the hexagonal symmetry and exhibits the linear energy dispersion. Finally, we 
explore the electronic properties of Si(lll)-surface structures on Ag surfaces. 

2. Calculations 

The total-energy electronic-structure calculations have been performed in the density 
functional theory using VASP code [36, 37]. As for the exchange-correlation functional, 
the vdW-DF [38, 39], being capable of treating the dispersion force is adopted [40]. 
Calculations for monolayer silieene and some of bilayer silieene have been also performed 
by the local density approximation (LDA) and the generalized gradient approximation 
(GGA) [13, 14]. The results obtained are essentially same among the three exchange- 
correlation functionals, although quantitative values for the total energies are different. 
The electron-ion interaction is described by the projector augmented wave method [41], 
and the cutoff energy of 250 eV in the plane-wave basis set is used. The Ag surface is 
simulated by a repeating slab model in which a five-atomic-layer slab is cleaved from 
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the face-centered-cubic (fee) Ag with the experimental lattice constant (4.09 A). The 
Ag slabs are separated from their images by the vacuum regions with the thickness 
of 14 A, 17 A, 21 A, and 24 A for the monolayer, bilayer, tri-layer and quad-layer Si 
deposited on Ag(lll), respectively. The geometry optimization is performed until the 
remaining forces become less than 0.02 eV/A. We have carefully examined the validity of 
the /c-point mesh in Brillouin zone (BZ) integration, and found that the spacing between 
adjacent fc-point being less than 0.016 A“^ is sufficient to obtain converged results. 

It has been recognized that the 3x3 and \/l x \/7 silieene are commensurate to 
the 4x4 and \/l3 x -^/Is Ag(lll), respectively [13]. To explore the 7r-bonded chain and 
DAS structures on Ag(lll) surface, we prepare initial structures by depositing the Si 
layers that are cleaved from the Si(lll) with the periodicity of the 2x2 (two times of 
the 2x1 unit cell) and 7 x 7 on the \/l x \/7 and 2\/^ x 2\/^ Ag(lll), respectively, 
where the lattice mismatches are 0.4 % and 1.4 %. Following the convention of the 
silieene layer, we define a single layer, i.e., a single bilayer, as consisting of 4 Si atoms 
at the upper positions and the remaining 4 Si atoms at the lower positions in the 2x2 
lateral cell of the Si(lll) surface. 

To assess relative stability among various stable and metastable structures, we 
introduce the cohesive energy Ec and the binding energy E]^. They are defined as 
A'c (A'Ag(iii) T Asi/isi F/tot)/Agi and E\) (A'Ag(iii) T F/gj—iayers A'tot)/ALsi, 
respectively. Here Atot, F'Ag(iii) and Esi-\ayeis are the total energies of the Si layers 
on Ag(lll), the clean Ag(lll) surface and the freestanding Si layers, respectively. Agi 
and Alsi are the total number of Si atoms and the number of Si atoms per Si layer, 
respectively, /igi is the chemical potential of Si which is adopted as the total energy of 
an isolated Si atom in this paper. The cohesive energy defined above is the energy gain 
to make Si ultra-thin layers on the Ag(lll) surface from the clean Ag surface plus the 
constituent Si atoms. This may be a measure of the catalytic ability of the Ag surface 
in forming the Si ultra-thin layers. The binding energy is the measure of the energy 
gain to put the freestanding Si structures on the Ag(lll) surface. In any case, E^ and 
Ab represent relative stability of various Si structures on the Ag(lll) surface. 

3. Stabilities and structures of Si layers on Ag(lll) 

In this section we present structural characteristics and energetics of ultra-thin Si atomic 
layers on Ag(lll) surfaces. We examine Si monolayer (ML), bilayer (BL), tri-layer (TL) 
and quad-layer (QL) structures. An interesting finding is the wealth of the structural 
multi-stability being composed of a variety of distinct structures with the cohesive- 
energy difference in the range of less than 100 meV per Si atom (meV/Si). The 7r-bonded 
chain and DAS structures observed on the Si(lll) surface are found to exist also in the 
ultra-thin Si layers on the Ag(lll) surface. 
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Figure 1. (color online) Top (upper panels) and side (under panels) views of geometry 
optimized structures for the 3x3 monolayer (ML) silieene on 4 x 4 Ag(lll). After 
optimization, the silieene shows the y/3 x \/3 (a) and the 3x3 (b) periodicities, 
respectively. The Si atoms in the top and bottom vertical positions are depicted by 
the large red and small blue balls, respectively. The large gray balls depict the positions 
of the substrate Ag atoms. The simulated lateral unit cell is indicated by the dashed 
(pink) lines in the top view, and its boundary is indicated by the solid (black) lines in 
the side views. The obtained rhombic y/3 x \/3 periodicity is indicated by the solid 
(orange) lines in (a). 


3.1. Monolayer Si on Ag(lll) 

Structures of ML silieene have been explored by our previous DFT calculations [13, 14]. 
In this subsection, we present those results and then add more details. We have found 
two distinct stable structures with 3x3 silieene periodicity, and four distinct stable 
structures with \/7 x \/7 periodicity. The calculated cohesive energies of all the six 
structures are in the range of 5.906 - 5.975 eV in LDA, inferring the possible observation 
of several distinct structures. Actually the calculated STM images of the three most 
stable structures agree with the experimental images [9, 10, 12] excellently and two of the 
three metastable structures have been recently observed in experiments [42]. The most 
stable 3x3 silieene obtained by the present vdW-DF calculation is shown in Fig. 1(b), 
where 6 of 18 Si atoms in a unit cell protrudes by 0.86 A, whereas the remaining 12 Si 
keep nearly the same height in the lower positions. The structural difference between 
the LDA and the vdW-DF is tiny. The cohesive energy obtained by vdW-DF is 5.410 eV 
(Table 1), which is between the LDA value (5.972 eV) and the GGA value (5.215 eV). 
We have also found a \/3 x y/S ML silieene structure [Fig. 1(a)] for the 3x3 sihcene/4 x 4 
Ag(lll) commensuration, in which one of six Si atoms in a unit cell protrudes about 1 
A, whereas the remaining hve Si keep nearly the same height in the bottom layer. Ghen 
et al. [16] have also found two distinct ^/3 x y/A structures both of which are different 
from the structure in Fig. 1(a). We have clarified that the structures obtained by Ghen 
et al. emerge when we perform the computation with insufficient k-point sampling [17]. 

The structural parameters for the stable x \/3 [Fig. 1(a)] and 3 x 3 [Fig. 1(b)] 
silieene on 4 x 4 Ag(lll) are shown in Table 1 along with the cohesive energies Ec and 
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Table 1. Calculated cohesive energy (eV/Si), binding energy E\y (eV/Si) and 
structural parameters for the silieene layers and the surface layers of Si(lll) on 
Ag(lll). dsi-Ag (A) is the spacing between the bottom Si layer and the topmost 
Ag layer. Note that the Si layer is buckled. Az (A) represents the averaged amount 
of the buckling of the top Si layer and dgi-si (A) is the interlayer distance defined as 
the minimum spacing between top Si layer and its underlying Si layer. 




Ec 


^si— Ag 

^si—si 

Az 

ML 

|co 

X 

5.387 

0.63 

2.37 

— 

1.05 


3x3 

5.410 

0.65 

2.17 

— 

0.86 

BL 

Vs X \/3 

5.339 

0.78 

2.25 

2.58 

0.96 


3 X 3-AA 

5.342 

0.79 

2.18 

2.47 

0.81 


V? X V7-AA' 

5.331 

0.80 

2.27 

2.39 

0.99 


2 X 1 

5.350 

0.58 

2.26 

2.25 

0.98 


7x7 

5.316 

0.63 

2.23 

1.40 

— 

TL 

\/3 X Vs 

5.384 

0.83 

2.26 

2.53 

1.08 


3 X 3-AAA 

5.364 

0.77 

2.18 

2.48 

0.89 


3 X 3-ABC' 

5.391 

0.85 

2.26 

2.41 

1.06 


V7 X V7-AAA' 

5.370 

0.80 

2.26 

2.40 

1.10 


2x1 

5.403 

0.81 

2.28 

2.24 

1.05 

QL 

Vs X Vs 

5.418 

0.86 

2.28 

2.52 

1.09 


3 X 3-AAAA 

5.384 

0.73 

2.18 

2.49 

0.88 


3 X 3-ABCA' 

5.423 

0.88 

2.26 

2.44 

1.17 


V7 X V7-ABCA 

5.404 

0.89 

2.26 

2.49 

0.95 


2 X 1 

5.434 

0.88 

2.27 

2.24 

1.06 


the binding energies E\y. From Table 1, the \/3 x \/3 structure is 23 meV/Si less stable 
than the 3x3 silieene on Ag(lll) which has been widely observed in experiments. This 
shows that the \/3 x \/3 ML silieene is hard to be synthesized, being consistent with the 
previous experimental results [19, 20]. As for the yr-bonded chain structure, commonly 
observed on Si(lll) surface, the topmost two surface layers show drastic reconstruction. 
Hence the ML silieene cannot evolve to the yr-bonded chain structure. 

3.2. Bilayer Si on Ag(lll) 

In this subsection, we present the results for the BL silieene. We focus on the 3x3 
BL silieene on 4 x 4 Ag(lll) surface and on the \/7 x \/7 BL silieene on \/l3 x ^/T3 
Ag(lll) , where the lattice mismatches are 0.4 % and 2.6 % (with respect to the (111) 
plane of diamond-structured Si), respectively. Another degree of freedom is the way of 
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Figure 2. (color online) Top (upper panels) and side (middle panels) views of the 
calculated stable structures and corresponding STM images (bottom panels) of the 
three rhombic \/3 x \/3 bilayer (BL) silieene on Ag(lll). The (a), (b) and (c) are 
labeled as strl, str2 and str3 of BL silieene, respectively. The large gray balls and the 
small yellow balls depict the positions of Ag atoms of substrate and the Si atoms of 
the bottom Si layer, respectively. The large red balls and small blue balls depict the 
positions of protruded and unprotruded Si atoms of the top Si layer. The lateral unit 
cells are indicated by the dashed (pink) lines in the top views, and their boundaries are 
depicted by the solid (black) lines in the side views. The obtained rhombic x \/3 
periodicity is indicated by the solid (orange) lines. 


stacking of the two Si layers. When the two planar hexagonal networks are stacked in 
the same way, we call it the AA stacking, whereas when the two networks are rotated 
by 60 degrees we call it the AB stacking, following the convention of graphite. From 
these two different stacking ways, we start the geometry optimization. We also examine 
stability of the vr-bonded chain structure and the DAS structure for this BL Si atomic 
layers on the Ag(lll) surface. One of the interesting features is the multi-stability of 
various structures of which the cohesive energies are in the range of 5.316-5.350 eV. 

Starting from the 3x3 BL silieene on the 4x4 Ag(lll) with AB stacking, we 
have reached three stable \/3 x ^/3 rhombic silieene structures shown in Fig.2. The 
three strnctures, labeled as strl, str2 and str3, are equivalent to each other when we 
neglect the existence of the Ag substrate: In the top layer, a single Si atom is protruded 
by about 1 A from the remaining five Si atoms in the lateral y/S x \/3 periodicity, 
similar to the -^/S x y/S structure of ML silieene in Fig. 1(a). In the bottom silieene 
layer, Si atoms are buckled with the amount of about 0.8 A, as they are on the (111) 
plane of diamond-structured Si. These features are common for the three structnres. 
Moreover, the total-energy difference among the three structures are within 0.36 meV/Si, 
inferring coexistence of all three structnres in usual experiments. These structnres can 
be transformed to each other by moving the rhombic \/3 x unit cells along their 
diagonal direction. 
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Figure 3. (color online) Top (upper panels) and side (middle panels) views of the 
calculated stable structures and corresponding STM images (bottom panels) of BL 
silieene and the 7r-bonded chain and the DAS structures on Ag(lll). (a), (b) The 
3x3 and the Cl x Cl BL silieene structures, respectively, (c), (d) The 2x1 and the 
7x7 BL Si(lll)-surface structures, respectively. The color codes are the same as in 
Fig. 2. The obtained 2x1 periodicity in (c) is indicated by the solid (orange) lines. 


The calculated STM images [43] of the three \/3 x ■sfz structures show the perfect 
rhombic geometry, which agrees well with that observed by Chen et al. in experiment 
below 40K [16]. As will be discussed in section 4, the flip-flop motion would happen 
between two of the three structures at higher temperatures. As a result, the honeycomb 
■\/3 X \/3 STM image [Fig. 8(a)] that is widely reported in the experiments [19, 21, 22] 
will be obtained. 

The structure of the 3x3 BL silieene on the 4x4 Ag(lll) with A A stacking keeps 
its periodicity even after the geometry optimization [Fig.3(a)]. The top layer of the 
structure resembles the 3x3 ML silieene on Ag(111) [Fig. 1(b)]. The corresponding STM 
image is further calculated , where each protruded Si atom in the top layer corresponds 
to a bright spot [bottom panel of Fig. 3(a)]. The STM image is nearly the same 
as that of the 3x3 ML silieene on Ag(lll) which has been extensively observed in 
experiments [9]. It is known that the usual experimental approaches (STM and LEED) 
can only reflect the geometry feature of the top silieene layer. This means that people 
are easily mislead by the STM/LEED geometries in identifying the thickness of silieene 
in experiments. Although, the height along particular direction of the silieene layers 
can be measured and thus used to distinguish silieene thickness, it does not work when 
the Ag surface is not well verified [44] or the Ag surface is not flat [20]. This indicates 
that not only the STM measurement but also the determination of the Si coverage is 
important for the structural identification. 

As for the s/l x Cl BL silieene on \/l3 x \/l3 Ag(lll), we have also considered both 
A A and AB stacking cases. The most stable structure we have found is the Cl x Cl 
BL silieene with AA' stacking [Eig. 3(b)], where A' means the in-plane dislodgment of 
the top Si layer from the proper A stacking. The calculated STM image [bottom panel 
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of Fig. 3(b)] further shows the perfect \/7 x \/7 periodicity with respect to that of 1 x 1 
silieene. The most stable AB-stacking and AA-stacking structures we have obtained 
present the 1x1 and \/7 x \/7 periodicities, which are 1.5 meV/Si and 29 meV/Si 
higher in energy than the AA'-stacking structure, respectively. 

The cohesive energies for the BL silieene on Ag(lll) are shown in Table 1, all 
of which are larger than that of freestanding silieene (4.75 eV/Si in our calculation), 
showing the catalytic function of Ag in forming silieene. On the other hand, the 
cohesive energies of BL silieene on Ag(lll) are smaller than those of the ML silieene 
on Ag(lll) by 37-94 meV/Si. This explains why the ML silieene structures are usually 
synthesized before the BL silieene in experiments [18, 19]. The 3x3 BL silieene with 
AA stacking on Ag(lll) [£'c=5.342 eV/Si, Fig. 3(a)] is the most stable among all the 
BL silieene structures, which is lower in the total energy by 3 meV/Si than the -s/S x y/3 
silieene structures (AB stacking. Fig. 2) widely observed in experiments [16, 18, 19, 20]. 
Considering that the 3x3 BL silieene exhibits nearly the same STM image as the 
3x3 ML silieene on Ag(lll) [9], this strongly indicates the existence of the 3 x 3 BL 
silieene structure on the Ag surface, which may be misrecognized to be the ML silieene 
in previous experiments. In addition, the x ^/7 BL silieene on Ag(lll) is higher in 
the total energy by about 10 meV/Si than the 3x3 and \/3 x \/3 structures, showing 
that it is rarely synthesized in experiment. 

To explore the possibility of Si(lll) reconstructed surface emerging on the Ag(lll) 
substrate, we deposit the topmost two surface layers of the vr-bonded chain structure [24] 
and the DAS structure [25, 26] on Ag(lll), though the two top-layer atoms in the DAS 
structure is incomplete for the full reconstruction in the DAS model. After geometry 
optimization, it is found that both structures are stable [Figs. 3(c) and 3(d)]. The 
obtained stable 2 x 2 Si surface shows the periodicity of 2 x 1 in which the bond network 
in the Si double layer contains five-membered and seven-membered rings and thus two 
top surface atoms in a unit cell form a 7r-bonded chain, being essentially identical to the 
TT-bonded structure on the Si(lll) surface. We have also calculated the corresponding 
STM image. As shown in the bottom panel of Fig. 3(c), the STM image agrees well 
with that of the pure 2x1 Si (111) surface [45]. These results clearly show the emergence 
of the 2x1 Si(lll)-surface structure on Ag(lll). 

The 7 X 7 Si(lll)-surface structure is also preserved when we deposit the top two 
Si layers on Ag(lll) [Fig. 3(d)]: the top layer composed of the 12 Si atoms in a unit 
cell keeps nearly the same geometry as that of the pure 7x7 Si(lll) surface, and the 
geometry of its underlying layer composed of the remaining 90 Si atoms is also basically 
preserved on Ag surface. The calculated STM image [bottom panel of Fig. 3(d)] is 
similar to that of the pure 7 x 7 Si (111) surface [46]. The above results clearly show 
that the bottom Si layer acts as a buffer layer on the Ag surfaces, which preserves the 
structure of the Si layers above it. 

We have also calculated the cohesive energies for the Si(lll)-BL system deposited 
on Ag(lll). It is found that the 7r-bonded chain structure on the Ag(lll) is the most 
stable among all the BL Si on Ag(lll). The cohesive energy of the yr-bonded chain 
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structure (£'c=5.35 eV/Si) is larger by about 10 meV/Si than those of the 3x3 (£'c=5.342 
eV/Si) and the -^/S x \/3 (£'c=5.339 eV/Si) BL silieene structures (Table 1). As we will 
explain below, this difference in the cohesive-energy mainly come from the difference in 
the stability of the top Si layer. 

A common featnre for these BL Si structnres (except for the 7x7 DAS structnre) 
on Ag(lll) is that half amonnt of Si atoms in the bottom layer located in the lower 
positions are chemically bonded with the Ag atoms, whereas the remaining half amount 
of Si atoms located at the higher positions are chemically bonded with the Si atoms 
in top layer (side views of Figs. 2 and 3). This makes the bottom layer of Si being 
sp^ bonded with the Ag surface atoms and the top-layer Si atoms. If we consider the 
Ag substrate and the bottom Si layer as a united system, it can be expected that the 
cohesive energies of this united systems are nearly the same for various BL Si structures 
on Ag(lll). Therefore, the cohesive-energy difference mainly comes from the snrface- 
energy difference of the top layer for all the BL Si structures on Ag(lll). To confirm 
this, we have explored the stable strnctures of pure Si(lll) surface with the \/3 x \/3 
periodicity and found that its cohesive energy is £'c=5.18 eV/Si, being smaller than 
the pnre 2x1 7r-bonded Si(lll) snrface by about 10 meV/Si. This amount coincides 
with the cohesive-energy difference between the \/3 x and 2x1 BL Si structures on 

Ag(lll). 

It shonld be mentioned that the 7x7 BL Si(lll)-surface strncture does not follow 
the rule stated above. The 7x7 DAS structure is most stable on Si(lll) surface, 
whereas it is least stable as the BL system on the Ag(lll) (Table 1). This is dne to 
the complex reconstruction in the DAS model in which subsurface atoms are relaxed 
substantially. The BL system on the Ag(lll) is not thick enough to incorporate the 
complex reconstruction. However, we expect that the 7x7 Si(lll)-surface structure on 
Ag(lll) would be the most stable when the number of Si layers increases. Yet this issue 
will not be discussed in this paper mainly due to the limited resource of the computation. 

3.3. Tri-layer Si on Ag(lll) 

In this snbsection, we present the resnlts for the TL silieene. Let us start with 3x3 TL 
silieene on 4 x 4 Ag(lll). Interestingly, we have fonnd three stable \/3 x \/3 rhombic 
silieene structures for the TL silieene [Fig. 4], as in the 3x3 BL silieene. The three 
strnctures are transformed to each other and the total energies of the three structures 
are nearly the same, showing the tristability of the ^/3 x y/S TL silieene. The three Si 
layers are stacked as ABC, being identical to the stacking of the Si(lll). The structnre 
of the top silieene layer of the TL silieene is very similar to that in the BL \/3 x \/3 
silieene. The calcnlated STM images are nearly the same with those of the \/3 x \/3 
BL silieene on Ag(lll) [the bottom panels of Fig. 4]. 

Since we have found the 3x3 AA-stacking BL silieene [Fig. 3(a)], we explore the 
similar 3x3 TL silieene strncture with AAA stacking where the top layer resembles that 
of the AA-stacking BL silieene. As a result, the 3x3 structure is found to be unstable and 
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Figure 4. (color online) Top (upper panels) and side (middle panels) views of the 
calculated stable structures and corresponding STM images (bottom panels) of the 
three rhombic x \/3 tri-layer (TL) silieene on Ag(lll). The (a), (b) and (c) are 
labeled as strl, str2 and str3 of TL silieene, respectively. The large gray balls and 
the small yellow balls depict the positions of Ag atoms of substrate and the Si atoms 
under the top Si layer, respectively. The large red balls and small blue balls depict the 
positions of protruded and unprotruded Si atoms of the top Si layer. The simulated 
lateral unit cell is indicated by the dashed (pink) lines in the top view, and its boundary 
is indicated by the solid (black) lines in the side views. The obtained rhombic \/3 x \/3 
periodicity is indicated by the solid (orange) lines. 



Figure 5. (color online) Top (upper panels) and side (middle panels) views of the 
calculated stable structures and corresponding STM images (bottom panels) of the TL 
silieene and Si(lll)-surface structures on Ag(lll). (a), (b) The two symmetry-related 
3x3 silieene structures with the AAA stacking, which are labeled as strl and str2, 
respectively, (c) The 3x3 silieene structure with ABC' stacking, (d) The a/7 x a/7 
silieene structure, (e) The 2 x 1 Si(lll) 7r-bonded chain structure. The color codes are 
the same as in Fig. 4. The obtained 2x1 periodicity in (e) is indicated by the solid 
(orange) lines. 
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relaxes into either of other two 3x3 structures with AAA stacking [strl, str2 shown in 
Figs. 5(a) and 5(b)], where 5, not 6 in the case of the BL silieene, of 18 Si toms in the top 
layer protrude by about 0.9 A. The resultant two structures are symmetrically equivalent 
if we neglect the existence of the Ag substrate. The total energy difference between the 
two symmetry-related structures is therefore nearly zero, showing the bistability of the 
3x3 structures in TL silieene. The calculated STM images of the two structures reflect 
the characteristics of the top-layer atom arrangements [bottom panels of Figs. 5(a) and 
5(b)]. 

We have found a new 3x3 TL silieene structure where the stacking of the three Si 
layers are dislodged [Fig. 5(c)]: The stacking is ABC' where the top layer C is dislodged 
in the lateral plane from the proper C stacking as in diamond-structured Si. In the top 
layer, 4 of 18 Si atoms are protruded by about 1 A, whereas the remaining 14 Si atoms 
locate in the lower positions. This causes very different STM profile from the other Si 
structures on Ag(lll) [bottom panel of Fig. 5(c)]. Since the cohesive energy of this 
structure is large enough (see below), this unusual STM pattern should be observed. 

As for the \/7 x \/7 TL silieene on ^/T3 x v^TS Ag(lll), we have found that the 
most stable structure keeps the \/7 x \/7 periodicity with the AAA' stacking [Fig. 5 
(d)]. The structure of top silieene layer and the STM pattern are nearly the same as 
those of the \/7 x BL silieene [Fig. 3 (b)]. 

The most prominent characteristic of the stable ML, BL, and TL silieene obtained 
in the present calculations is the protrusion of some of top-layer Si atoms causing the 
nearby buckling of the top layer. This structural relaxation was originally proposed to 
explain 2x1 Si(lll) surface by Haneman [47] before the yr-bonded chain model was 
proposed by Pandey [24]. Hence it is recognized that the restriction in the lateral plane 
caused by the Ag substrate realizes the buckling relaxation which is irrelevant in the 
free Si(lll ) surface. The buckling structure caused by the external environment may 
induce unexpected phenomena such as spin polarization [48]. 

To verify the stability of the vr-bonded structure in the TL silieene on the Ag(lll), 
we deposit the topmost three surface layers of the yr-bonded chain structure on Ag(lll). 
After geometry optimization, it is found that the yr-bonded chain structure is stable 
[Figs. 5(e)]. The calculated STM image [bottom panel of Fig. 5(e)] agrees well with 
that of the 2x1 yr-bonded chain structure on the Si(lll) surface [45]. 

The calculated cohesive energies for the TL Si structures on Ag(lll) are shown 
in Table 1. A common feature is that the cohesive energies in the TL structures are 
obviously larger than those of BL structures by tens of meV/Si, indicating that it is 
easier to obtain the TL Si structures on Ag(lll). The most stable silieene structure is 
the 3x3 silieene with ABC' stacking [£'c=5.391 eV/Si, Fig. 5(c)]. This is different from 
that in the BL silieene, where the 3x3 silieene with AA stacking is the most stable. 
The cohesive energy of the \/3 x y/S TL silieene structure (£'c=5.384 eV/Si) is very close 
to that of the most stable TL silieene structure, indicative of their coexistence. On the 
other hand, the 3x3 and \/7 x ^/7 silieene structures with AAA and AAA' stacking 
are smaller in Ec by 20-26 meV/Si than the most stable one, indicating that it may be 
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Figure 6. (color online) Top (upper panels) and side (middle panels) views of the 
calculated stable structures and corresponding STM images (bottom panels) of the 
three rhombic \/3 x \/3 quad-layer (QL) silieene on Ag(lll) labeled as strl (a), str2 
(b) and str3 (c), respectively. The large gray balls and the small yellow balls depict 
the positions of Ag atoms of the substrate and the Si atoms under the top Si layer, 
respectively. The large red balls and small blue balls depict the positions of the 
protruded and unprotruded Si atoms of the top Si layer. The simulated lateral unit cell 
is indicated by the dashed (pink) lines in the top view, and its boundary is indicated 
by the solid (black) lines in the side views. The obtained rhombic \/3 x a/ 3 periodicity 
is indicated by the solid (orange) lines. 


difficult to observe them in experiments. 

It is noteworthy that the 2x1 7r-bonded chain structure is stable even on the 
Ag(lll) substrate and further its cohesive energy (£'c=5.403 eV/Si) is larger than that of 
the most stable silieene structure (£'c=5.391 eV/Si) by about 10 meV/Si. This confirms 
our argument stated above: The relaxed or reconstructed Si(lll)-surface structures 
on Ag(lll) are more stable than the silieene structures owning to their lower surface 
energies. 

5^./. Quad-layer Si on Ag(lll) 

In this subsection, we present the results for the QL silieene. Starting from the 3x3 
TL silieene on the 4x4 Ag(lll) surface, we have reached the three stable a/S x \/3 
rhombic silieene structures, as in the BL and TL cases. The stacking is ABCA [Fig. 6] 
which seems to be a natural extension of the TL silieene with the ABC stacking. The 
structure of the top silieene layer and the STM pattern are nearly the same with those 
of the BL [Fig. 2] and TL [Fig. 4] \/3 x a/3 silieene on Ag(lll). 

The 3x3 silieene structures of AAAA [Fig. 7(a)] and ABCA' [Fig. 7(b)] stackings 
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Figure 7. (color online) Top (upper panels) and side (middle panels) views of the 
calculated stable structures and corresponding STM images (bottom panels) of the 
QL silieene and the 7r-bonded chain structures on Ag(lll). (a), (b) The 3x3 silieene 
structures with AAAA and ABCA' stacking, respectively, (c) The \/7 x Cl silieene 
structure with ABCA stacking, (d) The 2x1 7r-bonded chain structure. The color 
codes are the same as in Fig. 6. The obtained 2x1 periodicity in (d) is indicated by 
the solid (orange) lines. 


are also stable. The AAAA structure seems to be an extension of the structnres of TL 
silieene with the AAA stacking [Figs. 5(a) and 5(b)]. While the top-layer structure and 
the STM pattern present a different rotational angle with respect to the Ag surface. 
Other structure mapped by the symmetry operations is not fonnd for the QL case. In 
the ABCA' structure, the bottom 3 Si layers present the perfect ABC stacking. The 
top Si layer (A' layer) presents a in-plane dislodgment from the proper A stacking, the 
geometry of which resembles the C' layer of ABC'-stacking structure in the TL case, 
but with a different rotational angle with respect to the Ag snrface. Thus the STM 
pattern [bottom pattern of Figs. 7(b)] presents some similarity with that of the ABC' 
TL silieene on Ag(lll) [bottom pattern of Figs. 5(c)]. 

The most stable strncture of the \/7 x \/7 silieene on \/l3 x \/l3 Ag(lll) also 
presents the \/7 x -s/l silieene periodicity [ABCA stacking, Fig. 7(c)] as that in the 
TL case [Fig. 5(d)], where 3 of 14 Si atoms in the top layer protruded by about 1 
A constitute the triangle geometry. On the other hand, the STM images of the two 

X Cl strnctnres are different from each other [bottom panels of Fig. 7(c) and 
Fig. 5(d)] due to the different positions of the protruded Si atoms. This shows the Si 
thickness dependence of the most stable silieene structnres for the Cl x \/7 silieene on 
CT3 X ^/T3 Ag(lll). 

To verify the stability of the 2x1 vr-bonded structure in the QL silieene on the 
Ag(lll), we deposit the topmost four surface layers of the vr-bonded chain structure 
on Ag(lll) and perform the geometry optimization. As a result, the vr-bonded chain 
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structure is also stable [Fig. 7(d)] on Ag(lll) and the calculated STM image [bottom 
panel of Fig. 7(d)] agrees well with that of the 2x1 vr-bonded chain structure on the 
Si(lll) surface [45]. 

The calculated cohesive energies for the QL Si structures on Ag(lll) are shown in 
Table 1. The cohesive energies in the QL structures are obviously larger than those of 
the BL and TL structures, showing the increment of cohesive energies with the number 
of layers for multilayer silieene. The most stable silieene structure is the 3x3 structure 
with ABCA' stacking [£"2=5.423 eV/Si, Fig. 7(b)]. The cohesive energy of the \/3 x \/3 
silieene structure [£c=5.418 eV/Si] is very close to it, showing their coexistence. The 
%/7 X \/7 and 3x3 silieene structures with ABCA and AAAA stackings are smaller in Ec 
by 19-39 meV/Si than the most stable one, implying that it is hard to synthesize them in 
experiments. The cohesive energy (£c=5.434 eV/Si) of 2 x 17r-bonded chain structure is 
larger than that of the most stable silieene structure also by about 10 meV/Si, confirming 
that the relaxed or reconstructed Si(lll)-surface structures are more stable than the 
silieene structures on Ag(lll). 

The reason why the 2x1 yr-bonded chain and/or 7 x 7-DAS structures have not 
been observed on Ag(lll) in experiments is attributed to the synthesizing method, 
where the structures are always produced via depositing the Si atoms on Ag(lll). This 
makes the Si structures on Ag(lll) grown layer by layer: the BL Si structures are 
grown based on the structures of ML silieene and the TL Si structures are on the BL 
Si structures etc. Since the structure of ML silieene is very different from that of the 
Si(lll) layer, it is not expected to obtain the 2x1 BL Si(lll)-surface structure on 
Ag(lll) by using the synthesizing method widely adopted to make silieene. It is similar 
for synthesizing thicker (TL, QL and so forth) Si structures on Ag(lll). On the other 
hand, it is well known that the 2x1 Si (111) structure is obtained by cleaving the Si(lll) 
surface at room temperature, and it further converts to the more stable 7x7 structure 
when heated above 400 °C. This shows that different experimental approaches produce 
the different Si structures. Therefore, although the 2x1 and 7 x 7 Si structures are 
more stable than the silieene structures on Ag(lll), they have never been observed in 
experiment because no one has performed such experiments as heating the multilayer 
silieene structures at high temperatures to convert the silieene structures to the 2x1 
or 7 X 7 Si(lll) structures. 

Table 1 also shows the structure details and the binding energies for all the Si 
structures from ML to QL on Ag(lll). The binding energies for all the structures 
(the energy gain in the deposition of Si layers on the Ag surface) are in the range of 
0.6-0.9 eV/Si, larger than the typical vdW interaction energy manifested in the cases 
of graphene on metal surfaces by an order of magnitude [49], showing the covalent 
nature of the Si-Ag interface interactions. Moreover, the spacings between the bottom 
Si layer and Ag surface are in the range of 2.17-2.37 A. Considering that the atomic 
radii of Si and Ag are 1.18 and 1.65 A, respectively, the small spacing values confirm 
the formation of the Si-Ag covalent bonds. This result agrees with the recent DFT 
results which studied various monolayer silieene structures on Ag(lll) [31, 50]. While 
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Figure 8. (color online) Calculated reaction energy profile among the three rhombic 
y/3 X ^/3 structures of BL (a), TL (b) and QL (c) silieene on Ag(lll), shown in Figs. 
2, 4 and 6, respectively, (d) The calculated reaction energy profile between the two 
rotation-symmetric 3x3 structures of TL silieene with AAA stacking on Ag(lll) [Figs. 
5(a) and 5(b)]. The energy in the vertical axis is defined as {E — Egtri)/'^^ where E 
and Estri are the total energy of the structure at reaction coordinate and the total 
energy of strl, respectively, n is the number of atoms in the top silieene layer. The 
transition barriers between strl and str2, str2 and str3, and strl and str3 are depicted 
by the black squares, blue dots and red triangles, respectively. The insets show the 
STM images obtained by superimposing the images of strl and str2 structures. The 
simulated 3x3 lateral unit cells in the top views are indicated by the dashed (pink) 
lines. 


it disagrees with the argument in Ref. [51] that the Si-Ag interaction is weak. On the 
other hand, the interlayer distances between the top Si layer and the underlying layer 
for all the structures are smaller than 2.6 A. This shows the covalent nature of Si-Si 
interlayer interactions. As discussed below, the covalent interactions of Si-Ag interface 
and Si-Si interlayer would significantly influence the electronic properties of silieene. 

4. Flip-flop motion in multilayer silieene on Ag(lll) 

In this section, we present the flip-flop motion of multilayer silieene on Ag surface. We 
focus on the \/3 x \/3 tristable BL, TL and QL silieene structures and the 3x3 bistable 
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TL silieene structures. The calculated transition barriers are indicative of the flip-flop 
motion among/between the tristable/bistable silieene structures at low temperature. 
The results agree well with the STM images observed in recent experiments. 

In use of the nudged elastic band (NEB) method [52, 53], we have calculated the 
transition barriers among the three \/3 x \/3 structures, strl, str2, str3 of each BL, 
TL and QL silieene on Ag(lll) [Figs. 8(a), 8(b) and 8(c)]. The calculated transition 
barriers for the BL silieene are in the range of 7-9 meV/Si, whereas they are in the 
range of 14-28 meV/Si for the TL and QL silieene on Ag(lll). Anyhow, these barriers 
are low enough to allow the structural transition to be thermoactivated at temperatures 
of dozens of Kelvin. The insets of Figs. 8(a), 8(b) and 8(c) show superpositions of 
two of the three STM images for the BL, TL, and QL silieene structures, respectively, 
which present the perfect honeycomb V3 x \/3 structures. These results agree well 
with the recent experimental observations [16, 19, 20, 21, 22, 35]. It is noted that 
two other models had been recently proposed to explain the ^/3 x y/3 bilayer silieene 
structure observed in experiments [54, 55]. However, both of them can not explain all 
the experimental observations at lower (below 40K, rhombic structure) [16] and higher 
temperatures (above 77K, honeycomb structure) [16, 19, 20, 21, 22, 35]: The model 
in Ref. [54] cannot explain the honeycomb STM structure and the model in Ref. [55] 
cannot explain its conversion from the honeycomb STM structure to the rhombic STM 
structure at lower temperature [16]. From this sense, our model agrees better with the 
experimental observations so far. From the previous section, it is recognized that the 
bottom Si layer acts as a buffer layer on the Ag surfaces, above which the structure 
of Si layers can be well preserved. Thus the tristable \/3 x silieene structures 
and the low transition barriers among them are also expected to exist in the thicker 
multilayer silieene on Ag(lll). This explains the experimental observations where the 
same honeycomb \/3 x \/3 STM images are observed independent of the thickness of 
silieene [16, 18, 19, 20, 21, 22, 35]. 

The flip-flop motion also happens in other silieene structures. Fig. 8 (d) shows the 
calculated transition barriers between the two bistable structures of the 3x3 TL silieene 
with AAA stacking [Figs. 5(a) and 5(b)]. The very small energy barrier (3.5 meV/Si) 
shows that the flip-flop motion can happen at very low temperature. The inset shows 
the superposition of the STM images of the two structures, which exhibits the similar 
honeycomb 3x3 STM image as that of the 3x3 ML and BL silieene [bottom panel of 
Fig. 3(a)] on Ag surfaces. The results confirm the characteristic of flip-flop motion in 
multilayer silieene on Ag(lll) at low temperature. 

5. Dirac states and surface-localized tt states for the multilayer Si on 
Ag(lll) 

Previously, we have studied the electronic properties of 3 x 3 and \/7 x \/7 ML silieene 
on Ag(lll), which are very different from that of the freestanding silieene [13, 14]. A 
common feature for all the silieene structures on Ag(lll) is that the Dirac tt and tt* states 
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Figure 9. (color online) Calculated energy bands of the 3x3 silieene on 4 x 4 Ag(lll), 
where the silieene exhibits the \/3 x \/3 structure. The energy bands of the ML, BL, 
TL, and QL silieene on Ag(lll) are shown in (a), (b), (c) and (d). The energy bands 
of the freestanding ML, BL, TL, and QL silieene that are peeled from the Ag surface 
are shown in (e), (f), (g), (h), respectively. The origin of the energy is set to be £’f- 
The states that have characters of the tt (tt*) for the ML and BL silieene on Ag(lll) 
are indicated by the blue circles on T point. The surface- tt states in TL and QL silieene 
are indicated by the blue squares on T point. The states of the TL silieene on Ag(lll) 
that have characters of the surface- tt along the T — K direction are marked by the 
sequence of black squares in (c), which present linear energy dispersion. 
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Figure 10. (color online) Contour plots (left panels) and isovalue surface at its value 
of 10% of the maximum value (right panels) of the squared Kohn-Sham (KS) orbitals 
marked by the squares in Figs. 9(c) and (d) of the ^/3 x \/3 TL (a) and QL (b) silieene 
on Ag(lll). Figures (c) and (d) show the corresponding contour plots and isovalue 
surface of the KS orbitals marked in Figs. 9(g) and 9(h) of the \/3 x \/3 TL and QL 
silieene peeled from the Ag surfaces. The gray and blue balls depict Ag and Si atoms, 
respectively. 
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near Ep in the freestanding silieene disappear due to the strong silieene-Ag interactions. 
In this section, we present the results for the '/S x y/3 silieene on Ag(lll) because they 
are the only well observed multilayer silieene structures. Since the -s/S x \/3 structure 
has stable phases from ML to QL, such study can help us to understand the evolution 
of electronic properties with the increasing Si thickness on Ag surfaces. Additionally, 
we explore the effect of Ag substrate on the electronic properties of Si(lll)-surface 
structures, i.e., the 2x1 yr-bonded chain structure. 

The energy bands are calculated in the periodicity of 3 x 3 silieene on 4 x 4Ag(lll). 
Figs. 9(a) and 9(b) show the calculated energy bands for the ML and BL silieene on 
Ag(lll), where the K point in the BZ of the 1x1 silieene is folded on the F point. For 
both ML and BL cases, we have found no apparent energy band near Ep that shows 
the linear energy dispersion peculiar to Dirac tt (tt*) state. Instead, the states with 
characters of tt (tt*) appear deep below Ep [indicated by blue circles in Figs. 9(a) and 
9(b)]. The energy bands of ML and BL silieene which we obtain by peeling from the Ag 
surface are further calculated. It is found that the states with characters of tt (tt*) shift 
upward near Ep from the deep positions in the valence bands due to the annihilation of 
the interaction of Ag substrate [indicated by blue circles in Figs. 9(e) and 9(f)]. This 
shows that the covalent interactions between silieene and Ag surface makes the tt (tt*) 
states being deep in the valence bands. We also find no linear energy band peculiar 
to the Dirac tt (tt*) state near Ep for the TL and QL silieene on Ag(lll) [Figs. 9(c) 
and 9(d)]. We have additionally calculated the energy bands of the TL and QL silieene 
peeled from the Ag surfaces [Figs. 9(g) and 9(h)]. As a result, no Dirac characteristic tt 
(tt*) linear energy band near Ep is found. This implies that the covalent Si-Si interlayer 
interactions also cause the absence of Dirac states. These results strongly indicate the 
absence of Dirac states for multilayer silieene on Ag(lll) due to the covalent interactions 
at the Si-Ag interface and in the Si layers (Table. 1). 

On the other hand, we have found a new state near Ep for the TL and QL silieene 
on Ag surfaces, as indicated by the blue squares on F point in Figs. 9(c) and 9(d), 
which are mainly composed of the tt orbitals located at the top silieene layer [Figs. 
10(a) and 10(b)]. Here we call the new state the surface- tt state. The distributions 
of the Kohn-Sham (KS) orbitals for the surface-yr state [right panels of Figs. 10(a) 
and 10(b)] are similar with those of the Dirac tt (yr*) orbitals in the freestanding ML 
silieene, both of which present the hexagonal symmetry. Thus we can also expect that 
the surface-yr state exhibits linear energy dispersion near Ep [1]. Through the detailed 
analysis of the KS orbitals along the F-A' direction in BZ for all the energy bands, we 
have indeed observed such linear energy dispersion nearby Ep for the surface-yr state, 
as indicated by the black squares in Figs. 9(c) for the TL silieene on Ag(lll). The 
Fermi velocity of the surface-yr state is much smaller than that of the Dirac states in 
the freestanding ML silieene because of the (\/3 times) larger distances between the two 
nearest yr orbitals which cause the weaker transfer energy. Moreover, the KS orbitals of 
the surface-yr states are similar for the TL and the QL silieene on Ag(lll) [Figs. 10(a) 
and 10 (b]. This implies that the surface-yr state always exists near Ep for the \/3 x y/S 
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Figure 11. (color online) Calculated energy bands (a) and isovalue surface at its 
value of 25% of the maximum value of squared KS orbitals (side view) of the tri-layer 
2 X 2 Si(lll) on \/l x V? Ag(lll), where the Si surface presents 2x1 7r-bonded 
chain structure [Fig. 5(e)]. The energy bands and isovalue surface of squared KS 
orbitals of the TL Si structure peeled from Ag surface are further shown in (b) and 
(d), respectively. The states that have characters of the isolated tt (tt*) are indicated 
the blue squares on T point in (a) and (b). The gray and blue balls depict Ag and Si 
atoms, respectively. 


multilayer silieene on Ag(lll), irrespective of the thickness of silieene. We expect that 
the surface-TT state has great potential applications in the field of spin electronics [48]. 

To understand the origin of the surface-yr state, we have further calculated the 
energy bands of the x \/3 TL and QL silieene that are peeled from the Ag surface. 
As a result, the surface-yr states near Ep are also obtained [indicated by the blue squares 
in Figs. 9(g) and 9(h)] , although it contains yr orbitals on the subsurface Si layer for the 
TL silieene [Fig. 10(c)]. It is noteworthy that the geometries of the surface-yr orbitals 
of the TL and QL silieene peeled from Ag(lll) [Figs. 10(c) and 10(d)] are similar with 
those of the TL and QL silieene on the Ag(lll) [Figs. 10(a) and 10(b)]. This shows 
that the surface-yr state is originated from the unique y/S x y/3 Si surface structure. 
Therefore, one can always expect the appearance of surface-yr state near Ep in the thick 
\/3 X y/S multilayer silieene structures. 

Since the Si(lll)-surface structures are also stable on Ag surfaces, it is important 
to explore the effects of Ag substrate on the electronic properties of Si(lll)-surface 
structures. It is well known that the 2x1 yr-bonded chain Si surface presents 
a semiconductor property, where both the valence band maximum (VBM) and the 
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conduction band minimum (CBM) come from the isolated tt (tt*) states of the vr-bonded 
atomic chain [56, 57]. Here we have calculated the energy bands of the TL 2 x 1 vr- 
bonded chain structure on Ag(lll), where the unit cell adopted in the calculation has the 
periodicity of 2 x 2 Si(lll) on x Ag(lll) [Fig. 5(e)]. The calculated energy bands 
are shown in Fig. 11(a), which present the metallic property. Through the analysis of 
calculated KS orbitals, we hnd the metallic bands are mainly from the sates containing 
both the Ag orbitals and Si orbitals, as well as the sp band of Ag substrate. The states 
with characters of isolated vr and vr* appear 0.2 eV below and 0.5 eV above Ep on F 
point, respectively, as indicated by the blue squares of Fig. 11(a). The corresponding 
KS orbitals are shown in Fig. 11(c), which clearly shows the isolated vr and vr* states 
located on the vr-bonded atomic chain. 

We have also calculated the energy bands of the vr-bonded chain structure that are 
peeled from the Ag surface [Fig. 11(b)]. The isolated vr and vr* states are found about 
0.5 eV below and above Ep on F point [indicated by the blue squares in Fig. 11(b)], 
respectively. The corresponding KS orbitals are further shown in Fig. 11(d), which 
resemble those of the 2x1 structure on Ag(lll), although there are some additional 
vr orbitals on the subsurface Si layer, as shown in left panel of Fig. 11(d). The energy 
levels of the isolated vr and vr* states of the peeled 2x1 structure are close to those 
on the Ag surface, showing the minor effect of the Ag-Si interface interactions on the 
isolated vr and vr* states of top Si layer. We have also obtained the similar results for the 
QL vr-bonded chain structure on Ag(lll). Combining with the results of the s/A x \/3 
silieene on Ag(lll), it can be concluded that the strong Ag-Si interface interactions do 
not destroy the surface-related vr states of the multilayer Si structures on Ag surfaces. 

6. Conclusion 

We have performed the density-functional calculations for the ultra-thin atomic layers 
of Si on Ag(lll) surfaces, with the thickness of Si layers varying from monolayer to 
quad-layer. We have found several stable structures for the silieene on Ag(lll) with 
periodicities of s/Z x \/3, 3x3, and \/7 x \/7 with respect to the 1x1 silieene. We have 
also found that the 2x1 vr-bonded chain and 7x7 DAS Si(lll)-surface structures are 
stable on Ag surfaces, showing the crossover of silieene and Si(lll)-surface structures 
on Ag(lll). From the calculated cohesive energies, several common features can be 
deduced: (1) ML silieene is more stable than the BL and TL silieene on Ag(lll) 
but becomes less stable than QL silieene; (2) Cohesive energies of multilayer silieene 
and Si(lll)-surface structures on Ag(lll) monotonically increase with their thickness 
increasing; (3) The Si(lll)-surface structures are more stable than the silieene structures 
on Ag(lll) due to the lower surface energies. Moreover, we have found the structural 
tristability and bistability for the \/3 x \/3 and 3x3 multilayer silieene on Ag(lll), 
respectively. The calculated transition barriers are indicative of the flip-flop motion 
among/between the tristable/bistable structures at low temperature. The flip-flop 
motion between two of the three s/Z x s/Z structures produces the honeycomb STM 
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structures observed in experiments. A common feature for the multilayer Si structures 
on Ag surfaces is the covalent interactions at the Si-Ag interfaces and in the Si layers, 
which cause the absence of Dirac states in silieene. Instead, we have found a new state 
with the characters of vr orbitals located on the top silieene layer (surface-vr state), 
which appears near Fermi level and presents the linear energy dispersion. We expect 
the surface-TT states have great potential applications in the spin-related nano-devices. 
Finally, we have explored the electronic properties of 2 x 1 vr-bonded chain Si structure 
on Ag(lll). Our results show the robust characteristic of the surface-related tt states 
for the multilayer Si structures on Ag surfaces. 
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